Medium polarization isotopic effects on nuclear binding energies 
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There exist several eflective interactions whose parameters are fitted to force mean field predic- 
tions to reproduce experimental findings of finite nuclei and calculated properties of infinite nuclear 
matter. Exploiting this tecnique one can give a good description of nuclear binding energies. We 
present evidence that further progress can be made by taking into account medium polarization 
effects associated with surface and pairing vibrations. 
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Mean field theory is one of the most useful approxima- 
tions in all of physics and chemistry. In it we replace the 
many-particle Schrodinger equation by the single parti- 
cle Schrodinger equation. What started out as a problem 
too complicated to solve for all but the smallest systems 
becomes quite manageable when one deals with one par- 
ticle at a time (see e.g. [l|). Mass formulae based on 
mean field theory can achieve good agreement with ex- 
perimental findings. When phenomenological corrections 
are added to take into account selected correlation ef- 
fects (e.g. rotational energy correction, a Wigner term 
to take into account the proton-neutron interaction in 
nuclei with N nearly equal to Z) one can fit the experi- 
mental masses with a rms deviation of about 0.6 MeV. 

In the present paper we study some of the consequences 
zero point fluctuations associated with surface and pair- 
ing vibrations have on nuclear masses (in connection with 
the first subject see ref j^)- 



I. DETAILS OF THE CALCULATION 

In what follows we concentrate on a set of 121 spherical 
even-even nuclei (ranging from A = 1Q to A = 214). 

We started from a mean field calculated in the HF 
approximation, according to the procedure followed by 
Goriely et al. |3|. In particular, we made use of the 
Skyrme-type MSk7 interaction, adding a phenomenolog- 
ical Wigner term in order to take into account neutron- 
proton correlations in nuclei having nearly equal number 
of neutrons and protons. The pairing channel was treated 
in the BCS approximation using a zero-range interaction 



V^a^r{r,f') = Vo5[r-I') 



(1) 



in which we allow the pairing strength parameter Vo to 
be different for neutrons and protons and with an energy 
cutoff corresponding to one major shell [E — 41 ■ A~^'^ 
MeV). HFBCS solutions were computed in a box of 16 
fm radius, making use of a radial mesh of 0.1 fm. 



We then focused on the effect of correlations on the 
binding energies. The contribution associated with sur- 
face vibrations has been calculated in the quasiparticle 
random phase approximation (QRPA), according to the 
expression 



E'qrpa = -(2A +1)J2 ^aH J2 \Yilin)\ 



(2) 
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where hiux (n) are the energies of the n-th QRPA phonons 
of multipolarity A, Y^^{n) are the backwardsgoing am- 
plitudes (associated with the corresponding vibrational 
modes of energy hux {nj) and the sums over k and i run on 
QRPA quasiparticle indices (on RPA particle and hole in- 
dices for nuclei which are not superfluid). We don't con- 
sider here other important renormalization effects, e.g. 
self-energy effects, which affect the value of the effective 
mass (m* = 1.05m in the case of the MSk7 interaction) 
and the position of single-particle levels and level densi- 
ties. In QRPA calculations, using the same interactions 
adopted in the HF-t-BCS mean field, 2+ and 3~ vibra- 
tions have been calculated using a space dimension rang- 
ing from 40 MeV for the heaviest isotopes to 64 MeV for 
the lightest ones, exhausting the energy-weighted sum 
rule. Only the most collective low-lying QRPA phonons 
have been considered for the computation of correlation 
energies. The selection of the most collective phonons 
is achieved taking into account only those QRPA states 
with a reduced transition probability greater than 2% of 
non-energy weighted sum rule, see also ref. U ; for what 
concerns the individualization of the low-lying states, the 
fragmentation of the QRPA response imposes a case by 
case analysis, but a simple rule came out: quadrupolar 
(octupolar) states up to 7 MeV (5 MeV) contributed to 
the calculation of correlation energies, except for some 
nuclei (namely, Ar, Ca, Ti, Zr) for which we took into 
account quadrupolar (octupolar) states up to 10 MeV (7 
MeV). 

Pairing vibrations (see |a])l9) were treated in RPA. 
While it is true that well developed pairing vibrational 
bands based on doubly closed nuclei have been observed 
(e.g. in the case of ^°^Pb, see p), one expects pairing 



collectivity to be important also for systems with few 
nucleons outside closed shells. However, the correlation 
energy associated with pairing vibrations is calculated 
only on shell closures, while the correlations associated 
with pairing in open shell nuclei have been taken into 
account by means of the BCS approach. In any case, 
the above is an approximation which, if needed, can 
be removed from the calculation scheme. Single parti- 
cle levels have been computed using a Woods-Saxon po- 
tential with standard parametrization ITJ. A separable 
pairing interaction with constant matrix elements has 
been used to compute monopole pairing vibration con- 
tributions. In all calculations we have kept the contri- 
bution of only the lowest (n = 1) pair addition and pair 
subtraction modes, in keeping with the fact that, as a 
rule, the n ^ 1 modes are much less collective. The 
values of the associated pairing strengths in a specific 
nucleus have been computed as in ref. ^ (see note [g]). 
For neutron pairing vibrations the pairing strength has 
been computed in double closed shell nuclei and the re- 
sulting value has been used for isotopes with the same 
magic number of neutrons {^^0,^°Ca,'^^Gei,^^^Sn,'^°^Ph 
for N=8,20,28,82,126 respectively), while for N=50 the 
value has been computed in ^°Zr. For isotopic chains 
with magic atomic number, the correlation energy due to 
proton pairing vibrations has been calculated only for se- 
lected isotopes (i6-22o,40Ca,48Ca,"''Sn,i28Sn,208pb for 
Z=8,20,50,82) and then interpolated for the rest of the 
chain. 



the plot) are examples of magic nuclei which are under- 
bound with respect to near isotopes, while for N=132 
some opposite cases are observed. 
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FIG. 1: Deviations of computed mean field ground state en- 
ergies from experimental values as a function of neutron num- 
ber. The straight horizontal line is the E]"-'" - Ej'''' = line. 
Nuclei above this line are underbound. Isotopic chains are 
connected by solid lines. The full dots highlight calcium and 
lead isotopes, commented in the text. Figures come in color 
in the online version. 



II. RESULTS 



Mean field 



Fig. n shows the deviations of computed mean field 
(MF) ground state energies E'^"''-'^ from experimental val- 
ues E'^'^P {0,[l3)- The r.m.s. deviation 



rmsd 
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we obtain for this set of 121 nuclei is 0.724 MeV. Note 
that the parameters of the Skyrme force used are slightly 
different from the MSk7 developed by Goriely et al. in 
Q. We started with the MSk7 parametrization and, 
with the aid of a linear refit (which will be discussed 
below), we fitted the Skyrme parameters to the experi- 
mental binding energies of this set of 121 isotopes. The 
resulting parametrization was used as a starting point of 
the present work [llj . 

The shell structure dependence of the results is evident 
from Fig. ^ At variance with the results reported in ref. 
|3 obtained making use of the SLy4 Skyrme parametriza- 
tion, magic number nuclei are not, in our calculations, 
regularly overbound with respect to their neighbours. In 
fact, '*°Ca,'*^Ca and ^"*Pb (indicated by the arrows in 



B. Correlation energies 

Surface and pairing correlation energies are plotted 
in Fig. 121 Low-lying surface vibrations show, as ex- 
pected, a marked collectivity for all the isotopes calcu- 
lated. These corrections tend to bind magic nuclei less 
than their neighbours (compare Fig. |21 for the neutron 
shell closures at iV = 20,28,50,82). On the other hand 
pairing vibrations contributions show the opposite be- 
haviour, providing magic nuclei with an important extra 
binding energy. The sum of the two contributions is dis- 
played in the lower panel of Fig. [3 The differences which 
can be found between the results presented in ref. |j] and 
those presented in this work arise from the inclusion in 
the previous paper of quadrupole and hexadecapole pair- 
ing vibrations and from the fixing of some computational 
errors. 

After adding surface and pairing correlations to mean 
field binding energies, a refit of the Skyrme parameters 
is needed in order to minimize the rms deviation (PJ . Af- 
ter a first refit of the parameters one should recompute 
all nuclei in the mean field approximation and also re- 
compute the correlations. Since a full fit of the Skyrme 
interaction to the experimental binding energies is too 
costly, we will follow here the procedure of ref. [l3| . 
readjusting perturbatively the parameters of the force 




FIG. 2: Correlation energies associated to surface vibrations 
(SV, upper panel), pairing vibrations (PV, middle panel) and 
the sum of them (SV+PV, lower panel). Isotopic chains are 
connected by solid lines. All correlation energies are negative. 
The full dots highlight calcium and lead isotopes, as in FiglTI 



are sufScient to perform the refit of the interaction. Since 
certain hnear combinations of the Skyrme parameters are 
not determined by nuclear masses, one should not include 
them in the fit. According to 13] only four combinations 
are well fixed by the binding energies. The refit pro- 
cedure is linear in these integrals fi and one can vary 
the Ci coefficients in order to minimize the rms deviation 
between theoretical and experimental binding energies. 
We tested the linearity of the procedure recomputing the 
binding energies with the new parameters and comparing 
these energies to those obtained exploiting the linearity, 
finding excellent agreement. 

The results obtained after performing the refitting pro- 
cedure are shown in Fig. 13 The rms deviation for 
the refit of mean field plus surface and pairing corre- 
lations is 0.714 (middle panel of Fig. |3Jl, very close to 
the mean field value (0.724 MeV). It can be seen that 
adding the negative correlations has increased the devia- 
tion of the calculated masses from the experimental data 
for a few specific isotopic chains, which are now system- 
atically overbound (cf. for example the oxygen isotopes, 
indicated by an arrow in the middle panel of Fig. |2Jl, 
notwithstanding the refitting procedure. 



|l4| . This implies that the new interaction will be, in the 
parameters space, rather similar to the original MF inter- 
action. Thus, QRPA correlations are expected to change 
little under the refitting process. This is the reason why 
we leave unchanged the correlation energy contributions 
coming from QRPA. 

In the case of surface vibrations, we have compared 
our results with those recently obtained by Bender et al. 
[2,[l3l- We started from their mean field binding energies 
computed with the SLy4 Skyrme interaction within the 
HF-I-BCS framework for the 121 isotopes studied in the 
present paper; these binding energies display a rms de- 
viation from the experimental values of 2.45 (1.54) MeV 
before (after) a refit procedure is performed (see below for 
a brief explanation or [13 for a more detailed one) . If one 
computes the correlation energies obtained projecting on 
good angular momentum (J = 0), as the authors of [2| 
did, and adds these corrections to the mean field binding 
energies one obtains (after refitting) a rms deviation of 
1.45 MeV. If instead one adds the surface correlation en- 
ergies showed in the middle panel of Fig. |21 one obtains 
a rms deviation of 1.46 MeV. 

In what follows we give a brief description of the per- 
turbative refit performed (see ref. 13] for a detailed dis- 
cussion). The mean field Skyrme energy functional is 
decomposed into a sum of integrals fi, each of which is 
weighted by some linear combination c,; of the Skyrme 
parameters 
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E = Y,C^I^■ 



(4) 
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These integrals fi and the mean field binding energies 
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FIG. 3: Deviations of binding energies for mean field MF 
(upper panel), mean field MF corrected by QRPA correla- 
tions after the refit to experimental binding energies (middle 
panel) and mean field MF plus correlation energy fluctuations 
5E (discussed below in the text) after the refit to experimen- 
tal binding energies. Full dots highlight, as in previous fig- 
ures, calcium and lead isotopic chains. Isotopic chains are 
connected by solid lines. 

To gain perspective concerning this point, it is worth 
stepping back and looking at what happens when the 
ground state correlation energies are added to the iso- 
topic chains studied, without carrying out the refitting 
procedure. It is seen that, adding the correlation energy 
contributions, the deviation between theory and exper- 
iment for some of the nuclear chains studied displays a 



flatter isotopic dependence. An index of this flattening 
is provided by the rms deviation calculated with respect 
to the average deviation of the single isotopic chain. In 
other words, for the k — th isotopic chain 
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ak 



(5) 



where ak is the average deviation 
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The values of ak for all the 21 isotopic chains studied 
are displayed in Table ^ It can be seen that in many 
cases (13 out of 21) ak decreases compared to the mean 
fleld value, mostly due to the effect of pairing vibrations. 
The effect of surface vibrations alone is to increase ak 
compared to the mean field value in 16 chains out of 21. 
The total value of ak for all the isotopic chains, defined 
by the expression 



total 



' EkNkal 
EkNk 



(7) 



decreases by about 6% from 0.532 MeV for the MF cal- 
culation to the corrected mean field (0.502 MeV). 

In keeping with this finding, we have added to mean 
field only the fiuctuations of the QRPA correlations. 
That is, we added to each isotope the difference SE be- 
tween its correlation energy and the average correlation 
energy of the isotopic chain it belongs to. This proce- 
dure clearly cancels out the information related to the 
absolute value of the computed correlation energies, but 
keeps the, arguably, more important information con- 
cerning the isotopic dependence of these energies. In this 
context, we can remark that while the absolute value of 
the correlation energies associated with collective surface 
vibrations depends on the choice of the energy of the 
phonons which contribute to these correlation energies, 
the isotopic dependence is much more stable. This is 
because the properties of giant resonance modes change 
little along a given isotopic chain. Consequently, adding 
the fluctuations SE to MF selects the most stable fraction 
of the contributions associated with ground state corre- 
lation energies. The results obtained after refitting are 
shown in the lower panel of Fig. 13 The rms deviation is 
now 0.691 MeV. 



C. Further improvements 

Observing Table Q] one can realize that for some iso- 
topic chains (for example calcium, tin and sulphur) ak 
mantains nearly the same value or increases after adding 
surface and pairing vibrational correlation energies. The 
results shown in the Table arose from the application of 





# 


CTfe (MF) 


Crfe(MF + SV) 


(Tfc(MF + SV + PV(0+)) 


O (16 - 22) 


4 


0.430 


0.524 


0.436 


Si (30 - 36) 


4 


0.603 


0.642 


0.939 


S (30-38) 


5 


0.688 


0.850 


0.803 


Ar (36 - 46) 


6 


0.871 


0.946 


0.639 


Ca (40 - 48) 


5 


0.959 


1.153 


0.941 


Ti (40 - 54) 


8 


0.741 


0.751 


0.621 


Fe (50 - 58) 


5 


0.358 


0.468 


0.350 


Se (78 - 86) 


5 


0.136 


0.414 


0.250 


Kr (80 - 90) 


6 


0.296 


0.547 


0.298 


Sr (80 - 92) 


7 


0.548 


0.693 


0.519 


Zr (86 - 94) 


5 


0.535 


0.662 


0.461 


Mo (88-98) 


6 


0.352 


0.429 


0.208 


Sn (110-132) 


12 


0.498 


0.551 


0.473 


Te (124 - 136) 


7 


0.211 


0.316 


0.258 


Xe (130- 140) 


6 


0.333 


0.278 


0.301 


Ba (132-140) 


5 


0.373 


0.250 


0.354 


Ce (134 - 144) 


6 


0.612 


0.306 


0.407 


Sm (140-148) 


5 


0.543 


0.263 


0.333 


Hg (200-206) 


4 


0.270 


0.244 


0.416 


Pb (204-212) 


5 


0.518 


0.541 


0.417 


Po (206-214) 


5 


0.359 


0.502 


0.416 


O'total 


121 


0.532 


0.591 


0.502 



TABLE I: Rms deviations a^ of calculated ground state en- 
ergy from experimental values are shown. These deviations 
are calculated with respect to the average deviation for each 
isotopic chain (see eqs. |S]and|HJ. The Table displays: the 
mass number range (1st column), the number of isotopes cal- 
culated (2nd), the mean field rmsd ak (3rd), the rmsd ob- 
tained adding the contribution of surface vibrations to the 
mean field ground state energy (4th), the rmsd obtained 
adding the contributions from pairing vibrations as well (5th). 
The last row shows ototai for all the isotopic chains (see eq. 
UJ. All numbers are in MeV. 



the prescription described in section Q] This prescription 
is the same for all nuclei we included in this study, and 
this is in the spirit of a systematic study of the effect of 
these correlations to nuclear masses. Nevertheless it is 
worth looking at what we can gain by studying certain 
chains in more detail, in order to learn something that 
could be useful to improve our general prescription. If we 
focus on calcium and lead isotopes, for example, we know 
that in these nuclei well developed quadrupole pairing vi- 
brations have been observed Q- This lead us to include 
in these isotopes, besides the monopole pairing contribu- 
tion, also the quadrupole one. The corresponding results 
are shown in Fig. 01 and are summarized in Table Ull The 
improvement is evident. In the middle panel of Fig. [S] 
we show the results obtained performing the refit on all 
121 nuclei, adding quadrupole pairing vibrations for cal- 
cium and lead; in the lower panel, we show the results 
obtained adding only the flucutations 5E. In the latter 
case, the rms deviation for all 121 isotopes is 0.667 MeV 
(that is, a reduction of the rmsd from MF of nearly 8%). 

Analyses like this one could be pursued in other regions 
of the chart of nuclides, and it may be possible to find 
out some general rules that one can include in a general 



prescription. 
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FIG. 4: The effect of adding surface and pairing correlation 
energies to mean field binding energies are plotted for Ca 
(left) and Pb (right) isotopic chains. Solid lines represents 
mean field calculation results (MF), shown for all the 121 iso- 
topes in Fig. Q Dashed lines are obtained adding the surface 
correlations of the upper panel of Fig. |5| to the mean field 
binding energies (MF+SV), dotted lines are obtained adding 
the monopole pairing correlations of the middle panel of Fig. 
Has weU (MF-|-SV-|-PV(0+)), and dot-dashed lines represent 
the results obtained adding quadrupole pairing contributions 
as weU (MF-|-SV-|-PV(0+,2+)). The dashed horizontal lines 
indicate the average values of the first and the last curves, in 
order to better judge the improvement gained going from MF 
to MF+SV+PV(0+,2+). 





0"fe(MF) 


CTfc (MF + SV) 


Crfe(MF + SV+PV(0+)) 


O"fe(MF + SV+PV(0 + 


2+)) 


Ca 


().!•.'-,!• 


1.153 


0.941 


0.308 


Pb 


0.518 


0.541 


0.417 


0.296 



TABLE II: Rms deviations ak for calcium and lead after the 
addition of the correlation energy contributions described in 
the text are shown. These deviations are calculated with re- 
spect to the average deviation for each isotopic chain (see eqs. 
|5] and I^Jl . The Table columns display the mean field rmsd 
(2nd column), the rmsd for mean field plus surface correla- 
tion (3rd), the rmsd obtained adding monopole pairing con- 
tributions as well (4th) and the rmsd when also quadrupole 
pairing vibrational contributions are included. All numbers 
are in MeV. 



D. Influence of the various contributions 

It is worth looking at what happens when we add to 
the mean field binding energy only some of the contri- 
butions to correlation energies: data are shown in Ta- 
ble lllll It is evident, if one looks at the rms deviations 
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FIG. 5: The same for Fig. but with the inclusion of 
quadrupole pairing vibrational correlations for calcium and 
lead isotopes. 



obtained adding the fluctuations SE of the correlation 
energies (last column of Table IIIIjI , that the octupolar 
surface correlation energy has a small influence on the 
results. This is observed in all the three cases shown in 
Table Iml 

Comparing the first and the third rows of this Table 
(and analogously the second and the fourth), one can 
realize that the quadrupole contribution to pairing vi- 
bration correlations has a large influence on the results. 
In particular, adding this contribution to all isotopes, 
rather than adding them only to calcium and lead iso- 
topes, leads to a worse rms deviation (compare the first 
and the fifth rows of Table IIII|1 . 

The values of aj. (see eq. ISjl for the various prescrip- 
tions are shown in Table Hvl It is seen that the results in 
Table III and IV follow the same trends, the best agree- 
ment with experiment being obtained in the case (e), 
namely including quadrupole pairing vibrations for cal- 
cium and lead isotopes. 



III. SUMMARY AND CONCLUSIONS 

We have performed a systematic analysis of the contri- 
bution of medium polarization effects associated with col- 
lective surface and pairing vibrations on the binding ener- 
gies of even-even spherical nuclei. We have started from a 
HF-I-BCS mean field calculation based on a Skyrme-type 
interaction, obtaining a rms deviation of 0.724 MeV with 
respect to experimental masses. Then we have added 
the correlation energy associated with low-lying surface 
and pairing vibrational modes (calculated in the QRPA 
and RPA approximation respectively) which typically in- 
creases the calculated binding energies by a few MeV. 



MF+ 



sv 


PV 


rmsd 


rmsd 


2+ 3~ 


0+ 2+ 


MF+SV+PV 


MF+5E 


X X 


X X 


0.935 


0.788 


X 


X X 


0.942 


0.817 


X X 


X 


0.714 


0.691 


X 


X 


0.702 


0.694 


X X 


X only Ca 
and Pb 


0.729 


0.667 


X 


X only Ca 
and Pb 


1.225 


0.677 



(a) 
(b) 
(c) 
(d) 

(e) 
(f) 



TABLE III: Rms deviations obtained after adding some of 
the contributions to correlation energies to the mean field 
ground state of the 121 isotopes. The Table columns display 
the contributions taken into account in each case (first four 
columns), the rmsd (obtained after the refit) for mean field 
plus these correlations (5th column) and the rmsd (obtained 
after the refit) for mean field plus the fluctuations 5E of these 
correlations (6th). All numbers are in MeV. The values should 
be compared to the mean field rmsd of 0.724 MeV. 









CTk(MF + 


sv + PV 


) 






ak(MF) 


{a) 


(b) 


(c) (d) 


(e) 


(,/•) 


o 


0.430 


0.436 


0.442 


0.436 


0.442 


0.436 


0.442 


Si 


0.603 


1.113 


1.141 


0.939 


0.973 


0.939 


0.973 


s 


0.688 


0.831 


0.779 


0.803 


0.768 


0.803 


0.768 


Ar 


0.871 


1.047 


1.066 


0.639 


0.658 


0.639 


0.658 


Ca 


0.959 


0.308 


0.212 


0.941 


0.818 


0.308 


0.212 


Ti 


0.741 


0.587 


0.460 


0.621 


0.492 


0.621 


0.492 


Fe 


0.358 


0.281 


0.281 


0.350 


0.350 


0.350 


0.350 


Se 


0.136 


0.703 


0.789 


0.250 


0.161 


0.250 


0.161 


Kr 


0.296 


0.649 


0.776 


0.298 


0.319 


0.298 


0.319 


Sr 


0.548 


0.572 


0.718 


0.519 


0.550 


0.519 


0.550 


Zr 


0.535 


0.673 


0.814 


0.461 


0.502 


0.461 


0.502 


Mo 


0.352 


0.332 


0.522 


0.208 


0.245 


0.208 


0.245 


Sn 


0.498 


0.621 


0.525 


0.473 


0.436 


0.473 


0.436 


Te 


0.211 


0.413 


0.436 


0.258 


0.235 


0.258 


0.235 


Xe 


0.333 


0.596 


0.685 


0.301 


0.383 


0.301 


0.383 


Ba 


0.373 


0.674 


0.770 


0.354 


0.455 


0.354 


0.455 


Ce 


0.612 


0.667 


0.746 


0.407 


0.485 


0.407 


0.485 


Sm 


0.543 


0.573 


0.655 


0.333 


0.413 


0.333 


0.413 


Hg 


0.270 


0.663 


0.667 


0.416 


0.422 


0.416 


0.422 


Pb 


0.518 


0.296 


0.370 


0.417 


0.498 


0.296 


0.514 


Po 


0.359 


0.359 


0.295 


0.416 


0.344 


0.416 


0.385 


O'total 


0.532 


0.623 


0.657 


0.502 0.497 


0.464 


0.472 



TABLE IV: Rms deviations at of calculated ground state en- 
ergy plus correlation energies from experimental values are 
shown. These deviations are calculated with respect to the 
average deviation for each isotopic chain (see eqs. |S]and|SJ. 
The results in the different columns are obtained including 
different correlations, according to the labels reported in Ta- 
ble III. The results reported in the fourth column (labeled 
(c)) have already been listed in the last column of Table I. 
All numbers are in MeV. 



Subsequently we performed a linear refit of the parame- 
ters of the interaction in order to obtain the best agreee- 
ment with the experiment. Although this refit procedure 
had some difficulties in minimizing the rms deviation of 
binding energies (some isotopic chains turned out to be 
systematically overbound), we obtained in most isotopic 



chains a flattening of the deviations of the binding en- 
ergies from the experiment. Adding the fluctuations of 
the correlation energies, losing the information related to 
their absolute value, the linear refit was able to find a new 
set of Skyrme parameters which gives a rms deviation of 
0.691 MeV (that is, a reduction of about 4.5% from the 
mean field treatment). 

In order to improve the results achieved, we tried to 
better study some critical chains. For example the cal- 
cium isotopes, whose deviations from experimental bind- 
ing energies turned out to be insensitive to the addi- 
tion of the correlation energies considered, and which 
contribute in an important way to the total rmsd. We 
have been able to improve the description of calcium and 
lead isotopes including the correlations associated with 
quadrupole pairing vibrations. This reduced the rms de- 
viation for all 121 isotopes to 0.667 MeV (corresponding 
to about 8% reduction with respect to the MP rmsd). 
The reduction obtained is not strong, but it is sufficient 
to justify the effort made in going beyond the mean field 
approach. The improvement gained by the specific study 
performed on Ca and Pb isotopes can suggest that study- 
ing certain chains in more detail can allow us to learn 
something that could be useful to build a general pre- 
scription. 

Concerning surface vibrations correlations, we also 
compared our results with what was obtained by Bender 
et al. in ref. |2| with the projection on good angular mo- 
mentum and obtained a very good agreement, pointing 
to the question of how much the details of the calculation 
count. 

Our results show that, for the set of isotopes stud- 
ied, the degrees of freedom we investigated can help us 
in improving the mass formula built on a self-consistent 
mean field. However, our investigation was limited by 
the adopted refitting procedure, which made a perturba- 
tive survey of the parameters space, tying the results ob- 
tained to the starting Skyrme parametrization. Further 
investigations should adopt a more general refitting pro- 
cedure, exploring a larger region of the parameter space. 
Doing such a costly search for a better minimum of the 
rms norm would also imply that at every step one should 
re-calculate the mean field binding energies and also the 
correlation energies, because a larger change of the pa- 
rameters implies that these correlation energies can show 
a no more negligible variation. 

We conclude emphasizing that the aim of this work is 
not to achieve a mass formula with a very low rms devi- 
ation from experimental findings, but to investigate the 
effects of some degrees of freedom of the nuclear system 
which have to be taken into account (namely, surface and 
pairing vibrations). 
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